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Abstract 

We discuss the ground state and the low-lying excitations of the spin-half Heisenberg antiferro- 
magnet on the two-dimensional square-kagome lattice. This magnetic system belongs to the class 
of highly frustrated spin systems with an infinite non-trivial degeneracy of the classical ground 
state as it is known also for the Heisenberg antiferromagnet on the kagome and on the star lattice. 
The quantum ground state of the spin-half system is a quantum paramagnet likely with a finite 
spin gap and with a large number of non-magnetic excitations within this gap. The magnetization 
versus field curve shows plateaux as well as a macroscopic magnetization jump to saturation due 
to independent localized magnon states. 



1 



I. INTRODUCTION 



The magnetic properties of low-dimensional antiferromagnetic quantum spin systems have 
been a subject of many theoretical studies in recent years. These studies are motivated by 
the recent progress in synthesizing quasi-two-dimensional magnetic materials which exhibit 
e Xciti „ tUm effects Qflflfflfl. 

A lot of activities in this area were focused on frustrated spin-half Heisenberg an- 
iiferromagnets like the J\-Ji antiferromagnet on the square lattice (see, e.g. Refs. 



and references therein) and on the cubic 1J, [l5( lattice, the Heisen- 



berg antiferromagnet (HAFM) on the star lattice 
on the kagome lattice (see the reviews 
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16 



171 ] and last but not least the HAFM 
] and references therein). Due to the 



extreme frustration the HAFM on the kagome and the star lattices shows an infinite non- 



trivial degeneracy of the classical ground state. Furthermore, both 



16 



spin lattices exhibit a 



13, 



21 



Although 



magnetization jump to saturation due to localized magnon states 
there is most likely no magnetic ground state order for the quantum spin-half HAFM on 
both lattices, the nature of both quantum ground states and also the low-lying spectrum are 



basically different. It was argued 



16| that the origin for this difference lies in the existence of 



non-equivalent nearest-neighbor (NN) bonds in the star lattice whereas all NN bonds in the 
tagome lattice are equivalent. Another striking difference relevant for magnetic properties 
23j] lies in the number of spins in the unit cell which is odd for the kagome lattice but even 

for the star lattice. As a result of the interplay between quantum fluctuations and strong 



frustration for the kagome lattice the quantum ground state is a quantum s pin 



iquid with 



very short-ranged spin, dimer, and chirality correlations (see e.g. [itI Il8l I20L l24l |2£ 
(small) spin gap to the triplet excitations and an exceptional density of low-lying singlets 
below the first magnetic excitation. On the other hand, for the star lattice one meets a 
so-called explicit valence-bond crystal with a well-pronounced gap to all excitations which 
can be attributed to the non-equivalence of the NN bonds and to the even number of s = 1/2 
spins in the unit cell 16J, Il2 1 • 



27 



lattice 



In this paper we consider the spin-half HAFM on the square-kagome [26 , 
(see Fig.[TJ). The square-kagome lattice is built by regular but also by non-regular polygons 
and it has two non-equivalent sites. Therefore it does not belong to the class of so-called 



uniform tilings 0,|2£| (like, e.g. square, triangular, star or kagome lattice). Nevertheless, 
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FIG. 1: The square-kagome lattice with N = 30 sites. The two topologically inequivalent nearest- 
neihghbor bonds are distinguished by solid (square bonds J5) and dashed lines (triangular bonds 
Jt). 

there exist some important geometrical similarities to the kagome and also to the star lat- 
tices. Similar to the kagome lattice it has coordination number z — 4, the even regular 
polygons (hexagons for the kagome, squares for the square-kagome lattice) are surrounded 
only by odd regular polygons (triangles) and both lattices contain corner sharing triangles. 
As a result the HAFM on the square-kagome lattice is also strongly frustrated and exhibits 
an infinite non-trivial degeneracy of the classical ground state. The similarity to the star 
lattice consists in the existence of non-equivalent NN bonds and in the fact that both lat- 
tices have an even number of spins in the unit cell. Moreover the classical ground state 
of the HAFM on the star lattice also exhibits an infinite non-trivial degeneracy. Due to 
these similarities we can expect that the HAFM on the square-kagome lattice is another 
candidate for a quantum paramagnetic ground state. However, the question arises, whether 
the quantum ground state displays similar properties as that for the kagome lattice or as 
that for the star lattice or none of them. 
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II. THE MODEL 



The geometric unit cell of the square-kagome lattice contains six sites and the underlying 
Bravais lattice is a square one (see Fig. |1J). For this lattice we consider the spin-half HAFM 
in a magnetic field h 

H = J ^ Si ■ Sj — hS z , (1) 
to) 

where the sum runs over pairs of neighboring sites (ij) and S z = ^ Sf. As mentioned above 
the square-kagome lattice carries topologically inequivalent NN bonds J5 (square bonds, 
solid lines in Fig. |TJ) and Jt (triangular bonds, dashed lines in Fig. [TJ see also Fig. |2J). For 
the uniform lattice these bonds are of equal strength Js = Jt = J and we set J = 1 in what 
follows. 



III. SEMI-CLASSICAL GROUND STATE 



In the classical ground state for h = the angle between neighboring spins is 2n/3. Since 
the triangles are "corner sharing", there is a non-trivial infinite degeneracy resulting from 
the possible rotation of two spins on a triangle (see also Fig. |2J). The classical ground state 
energy per bond is eg lass = —0.125 assuming classical spin vectors of length s = 1/2. Similar 
to the kagome and the star lattices there are two variants of the classical ground state, 
shown in Fig. El being candidates for possible magnetic ground state ordering. 

To discuss the influence of quantum fluctuations on a semiclassical level we perform a 
linear spin-wave theory (LSWT) starting from the coplanar classical ground states. We 
have to consider six types of magnons according to the six sites per unit cell. As for the 
kagome 0, [ill Q and the star lattice 16j the spin-wave spectra are equivalent for all 
coplanar configurations satisfying the classical ground state constraint. We obtain six spin- 
wave branches, three optical branches, one acoustical and two dispersionless zero modes. 
Thus also flat zero modes appear as it is observed for the kagome and star lattice case. 
There is no 'order-by- disorder' selection among the coplanar classical ground states due to 
the equivalence of the spin-wave branches obtained from LSWT, exactly like for the kagome 
lattice 0, [yj and the star lattice 

The ground state energy per bond for s = 1/2 in the LSWT is Cq = —0.236555. Due to 



the flat zero modes the integral for the sublattice magnetization diverges 



32j which might be 



is. 
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FIG. 2: Two variants of the ground state of the classical HAFM on the square-kagome lattice: 
The state on the left side has a magnetic unit cell which is three times as large as the geometric 
one and resembles the v3xv3 state of the kagome and the star lattices. For the state on the 
right side the magnetic unit cell is identical to the geometric one and corresponds to the q=0 state 
of the kagome and the star lattices. The dotted ellipses show further degrees of freedom of the 
highly degenerate classical ground state. 

understood as some hint for the absence of the classical order. Although on the semiclassical 
LSWT level the square-kagome, the kagome and the star lattices exhibit almost identical 
properties, the situation will be changed taking into account the quantum fluctuations more 
properly. 



IV. EXACT DIAGONALIZATION 



To take into account the quantum fluctuations going beyond the semiclassical LSWT 
we use Lanczos exact diagonalization (ED) to calculate the ground state and the lowest 
excitations for the s = 1/2 HAFM at h = on finite lattices of iV = 12, 18, 24, 30, 36 sites 
with periodic boundary conditions. The ground states of all those systems are singlets and 
the ground state energy per bond e and the degeneracy of the quantum ground state das 
are given in Table UJ Furthermore, we give in Table H] the gap to the first triplet excitation 
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(spin gap) A. Note that eo and A are significantly smaller than the corresponding values 
for the star lattice but of comparable size as the values for the kagome lattice. 

Now we compare the spin-spin correlations with those for the HAFM on the triangular, 
kagome and star lattices in FigEl For the triangular, kagome and star lattices we consider 
the strongest correlations as a measure for magnetic order for the largest finite lattices 
accessible for ED and present in Fig. [21 the maximal absolute correlations | (SfSj) | max for 
a certain separation R = |R$ — Rj| versus R. Contrary to those lattices the square- 
kagome lattice contains two inequivalent sites. Hence we present for the square-kagome 
lattice all different correlations \(SfSj)\ in Fig El Note further that we prefer to present 
the correlations for the finite square-kagome lattice with N = 30 sites, since it has better 
geometrical properties than the largest square-kagome lattice considered (N = 36). As 
expected we have very rapidly decaying correlations for the disordered kagome and star 
case, whereas the correlations for the Neel ordered triangular lattice are much stronger for 
larger distances and show a kind of saturation for larger R. The decay of the correlations 
for the square-kagome lattice is also very rapid thus indicating the lack of long-range order 
in the spin-spin correlation function. The two non-equivalent NN bonds carry very similar 
spin correlations, its difference for iV = 30 is only about 10%, which is in contrast to the 
star lattice where the two non-equivalent NN bonds differ by a factor of 3.5 [16] . 

Let us now discuss the low-lying spectrum of the star lattice (see Fig. Hj), following the 
lines of the discussion of the spectrum for the triangular j^j], the kagome lattice 24, |^ 



n 

and the star lattice [lfij . The lowest states E^^S) shown in Fig. |U are not well described 
by the effective low-energy Hamiltonian H e g ~ Eq + S 2 /2Nxo of a semiclassically ordered 
system: One can see rather clearly that the dependence E min (S) vs. S(S + 1) is not a 
linear one and there are no separated so-called quasi degenerate joint states 3^ which in 
the thermodynamic limit could collapse to a ground state breaking the rotational symmetry. 
Note further that the symmetries of the lowest states in each sector of S cannot be attributed 
to the classical ordered ground states shown in Fig. El These features are similar to the 
ha g „ m e QQ and the star Q lattices. Bet there is one stn k in g difference between the 
kagome lattice and the star lattice. While the former one has an exponentially increasing 



number non-magnetic singlets filling the sin. 
singlets were found for the star lattice jli] . \24 



et-triplet gap (spin gap) no such low-lying 



251 ] . This difference was attributed to the 



non-equivalence of NN bonds in the star lattice and the resulting dimerization of the ground 
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FIG. 3: The absolute value the spin-spin correlations |(SjSj)| versus R = |R$ — Rj| for the HAFM 
on the square- kagome (N = 30), the kagome (JV = 36), the star (JV = 42) and the triangular 
(N = 36) lattices. For the kagome, star and triangular lattices we present only the maximal values 
of |(SjSj)| for a certain separation R (the lines are guides for the eyes), for the square-kagome 
lattice we present all different values for |(SjS 3 -)| obtained by averaging over the four degenerate 



ground states. Note that the data on for the kagome lattice coincide with those of Ref. |3Jj| and 



the date for the triangular and the star lattice with those of Ref. 



y 



.y 



state. Though the square-kagome lattice has also non-equivalent NN bonds its spectrum 
is different from that of the star lattice, rather it shows similar to the kagome lattice a 
large number N s of non-magnetic excitations within the singlet-triplet gap. We find N s = 6 
(JV = 12), 13 (JV = 18), 17 (JV = 24), 47 (JV = 30), 38 (JV = 36). These numbers increase 
with growing size (except for JV = 36, which might be attributed to the lower symmetry 

1, where an 



of this finite lattice) but are smaller than those for the kagome lattice 
exponential increase of N s with JV was suggested. Our data for the square-kagome lattice 
do not allow a secure conclusion about a possible exponential increasing of N s with JV. 
For the discussion of magnetic long-range order we use the following finite-system order 
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FIG. 4: Low-energy spectrum for the HAFM on the square- kagome lattice (N = 30) (the inset 
shows the k points in the Brillouin zone). 



parameter 



17] 



in • = ( — > ! , (2) 

which is independent on any assumption on eventual classical order. The value m^ ass for the 
two ordered classical ground states shown in Fig. El is m^ ass = \ \/2/3, which is the same as 
for the classical v^xv^ and q=0 states on the kagome and on the star lattices. 

The numerical values for (m + ) 2 are collected in Table |U The values of (m + ) 2 for the 
square-kagome lattice are comparable to those for the kagome lattice but are slightly smaller 
than the corresponding values for the star lattice jlfi| . 

To estimate the values of eo, A and m + for the infinite square-kagome lattice we have 
extrapolated the data from Table H] to the thermodynamic limit according to the standard 



formulas for the two-dimensional spin-half HAFM (see, e.g. [17, 35, 36]), namely eo(iV) 
e (oo) + A 3 N~^ + 0(N~ 2 ) for the ground state energy per bond, m + (N) = m + (oo) + 
B X N^ + OiN- 1 ) for the order parameter, and A(N) = A(oo) + G 2 N~ 1 + 0{N~^) for the 
spin gap. In Table ITT1 the results of these extrapolations are presented and compared to those 
obtained for spin-half HAFM on the kagome and on the star lattices. Our data suggest a 



8 



TABLE I: Ground state energy per bond ep, ground state degeneracy dcs, spin gap A and square 
of the order parameter (m + ) 2 of the spin-half HAFM on finite square- kagome lattices. 

N 12 18 24 30 36 

eo (das) -0.226870 (1) -0.223767 (2) -0.224165 (1) -0.221527 (4) -0.222197 (3) 
A 0.382668 0.290191 0.263906 0.188865 0.139550 

(m + ) 2 0.184160 0.116455 0.086735 0.068618 0.060475 



TABLE II: Results of the finite-size extrapolation of the ground state energy per bond eo, the 
order parameter m + and the spin gap A of the spin-half HAFM on the square- kagome lattice. For 



comparison we also show results for the kagome and star lattices taken from Ref. 



16 



17. To see the 



effect of quantum fluctuations we present m + scaled by its classical value w-^ ass for the two ordered 
states shown in Fig. [21 (The negative, but very small, extrapolated values for the square-kagome 
and the kagome lattices are an artefact of the limited accuracy of the extrapolation. We interpret 
these negative values as vanishing order parameters.) 



lattice square-kagome kagome star 
e -0.2209 -0.2172 -0.3091 

A 0.052 0.040 0.380 

m+/m+ ass -0.032 -0.036 0.122 



small but finite spin gap and a vanishing order parameter. 

The values of the extrapolated quantities of the square-kagome lattice are very close to 
those of the kagome lattice. Therefore these data clearly yield evidence for a magnetically 
disordered quantum paramagnetic ground state of the spin-half HAFM on the square- 
kagome lattice which is most likely similar to that of the kagome lattice. 
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FIG. 5: Magnetization curves of some finite spin-half HAFM systems on square-kagome lattice 
in a magnetic field (N = 24, 30, 48, 54). 

V. MAGNETIZATION PROCESS 



In this Section we briefly discuss the magnetization versus field curve for some finite 
square-kagome lattices. The magnetization m is defined as m = 2(S Z )/N. We focus on 
those finite lattices having optimal lattice symmetries, i.e. N = 24, 30. In the high field 
sector we are able to present also data for N = 48 and iV = 54. 

The results are shown in Fig. Due to the spin gap (see Tables |l] and |H} one observes 
a small zero- field plateau. Clear evidence for a further plateau is found at m = 1/3 which 
can be attributed to the presence of tr iang les j^]. Note that a m = l/3_plateau is also 



observed for the triangular 



37 



38l . |39(, the kagome 



17 



21 



37 



40 



dp J 





and the star 



171 ] lattices. 



At the saturation field h s — 3 a jump 
this jump was discussed already in Refs. 



in the magnetization curve appears. The presence of 



28, 



42j and is related to the existence of independent 
localized magnon states found for a class of strongly frustrated spin lattices ^ 21, 13 



among them the kagome and the star lattices. In the case of the square-kagome lattice 
these localized magnons live on the squares. The height of the jump is 5m and is related to 
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the maximum number n max of independent localized magnons which can occupy the lattice. 
For the square- kagome lattice we have n max = iV/6 and consequently 5m = 1/3. We 
mention, that these localized magnon states are highly degenerate thus leading to a finite 



ry q< 

H 



residual T = entropy at the saturation field h s = 3 [13, |43j. Just below the jump, i.e. at 
m = 2/3 there is evidence for another plateau. Its width was estimated in Ref. |28| by finite 
size extrapolation to Ah »s 0.33 J for the infinite system. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have discussed the ground state properties of the spin-half Heisenberg 
antiferromagnet on the square-kagome lattice. This lattice has similarities with the 
kagome as well as with the star lattice. The kagome and the square-kagome lattices have 
coordination number z = 4 and are built by corner sharing triangles. The star lattice 
(z = 3) shares with the square-kagome lattice the property to have two non-equivalent 
nearest-neighbor bonds and to have an even number (namely six) sites per unit cell (note 
that the kagome lattice has three sites per unit cell and all nearest-neighbor bonds are 
equivalent). On the classical and on the semiclassical level of linear spin wave theory 
the ground state of the Heisenberg antiferromagnet on all three lattices exhibits very 
similar properties. However, it was argued that in the extreme quantum limit s = 1/2 
just these geometrical properties of the star lattice in common with the square-kagome 
lattice but different to the kagome lattice lead to different quantum ground states for 
the star and the kagome lattices. Interestingly, our results for the square-kagome 
lattice lead to the conclusion that the quantum ground state of the Heisenberg antifer- 
romagnet on the square-kagome lattice is similar to that of the kagome lattice. We 
find evidence for a spin-liquid like ground state with a small gap of about J/20 and 
a considerable number of low-lying singlets within this spin gap. Contrary to the star 
lattice case we do not see here a tendency towards forming a valence bond crystal ground 
state. The magnetization curve of the s = 1/2 HAFM on the square-kagome lattice 
shows a jump just below saturation and three plateaux at m — 0, 1/3 and 2/3. The 
low-energy excitations present near saturation field promise large magnetocaloric effects . 
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